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Abstract: In fluid dynamical models the freeze out of particles across a three dimensional space-time hypersurface 
is discussed. The calculation of final momentum distribution of emitted particles is described for freeze out surfaces, 
with both space-like and time-like normals, taking into account conservation laws across the freeze out discontinuity. 


1. Introduction 

The freeze out of particle distributions is an essential part of continuum or fluid dynamical reaction models. From 
the point of view of observable consequences this is one of the most essential parts of the model. On the other hand 
this step is not based on fluid dynamical principles and governed by a large variaty of ad hoc assumptions. The 
freeze out can be considered as a discontinuity across a hypersurface in space-time. 

The general theory of discontinuities in relativistic flow was not worked out for a long time, and the 1948 work 
of A. Tautu discussed discontinuities across propagating hypersurfaces only (which have a space-like normal vector, 
da^da^ = —1). Events happening on a propagating, (2 dimensional) surface belong to this category. 

Another type of change in a continuum is an overall sudden change in a finite volume. This is represented by a 
hypersurface with a time-like normal, da^da^ = +1, called confusingly both space-like and time-like surface in the 
literature. In 1987 Taubs approach was generalized to both types of surfaces,□ making it possible to take into account 
conservation laws exactly across any surface of discontinuity in relativistic flow. This approach also eliminates the 
imaginary particle currents arising from the equation of the Rayleigh line. When the EoS is different on the two sides 
of the freeze out front these conservation laws yield changing temperature, density, flow velocity across the front. 

In fact the freeze out surface is an idealization of a layer of finite thickness where the frozen out particles are 
formed, and the interactions in the matter become gradually negligible. The dynamics of this layer can be described 
in different kinetic models or four-volume emission models .0 The zero thickness limit of such a layer is the idealized 
freeze out surface. 

The invariant number of conserved particles (world lines) crossing a surface element, da^ , is dN = da^ , and 

the total number of all the particles crossing the FO hyper-surface, S, is N = j s da M . This total number, N, 
and the total energy and momentum are of course the same at both sides of the freeze out surface. If we insert the 
kinetic definition of iV M 

N>i = J P>1 fFo(x,p-,T,n,u v ) , 

into the integral over the freeze out hypersurface, S, we obtain the Cooper-Frye formula:!!! 

E'cfip = J f FO ( x ’ p ’ T ’ n ’ u ^ > (!) 

where /po(:r,p; T, n, u u ) is the post FO phase space distribution of frozen-out particles which is not known from 
the fluid dynamical model. Problems usually arise from the bad choice of this distribution. First of all, to evaluate 
measurables we have to use the correct parameters of the matter after the FO discontinuity! 

If we know the pre freeze out baryon current and energy-momentum tensor, Nq and Tq V , we can calculate localUj, 
across a surface element of normal vector da 11 the post freeze out quantities, N 11 and , from the relations!™: 
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[N^ da M ] = 0 and [T^" da^] = 0, where [A] = A — A 0 . Iil numerical calculations the local freeze out surface can 
be determined most accurately via self-consistent iterationE’EI This fixes the parameters of our post FO momentum 
distribution, fFo{x,p-,T,n,u v ). 

For example we can illustrate the effect of conservation laws for a situation where the frozen out matter is massless, 
baryonfree Bose gas. Then, the conservation laws across the freeze-out surface with timelike normal vector da^ are 
[T^da^ = 0, In the most general (three dimensional) case there are four parameters to be determined from the 
conservation laws: The final, post FO temperature, T, and three components of the velocity, u. The energy- 
momentum tensor on the pre freeze-out side, and the normal to the surface are given. The post freeze-out energy- 
momentum tensor is of the form T^ v = (e -\-p)u^u 1 ' — pg ,lv , where the energy density, pressure, and temperature are 
connected by the EoS: e = asBT 4 = 3 p, where asB is the Stefan-Boltzmann constant. Then T^ v = (e+p)w /i u I/ — pg^, 
can be written as a vector equation: 

{Au^u v da v — da = xa M , (2) 

where 


x = = Tq V da v . 

Taking the normal projection of (§ and the norm of the four-velocity, the solution for the four quantities we are 
looking for will bei given by: 


a M dcr M + {a^da^) 2 + 3 a^a M 


xa M + da M 
2 ^Jxa^da^ + 1 


( 3 ) 


a. Idealized freeze out across propagating discontinuities. One can go a step further in the study of freeze out 
process. We usually assume that the pre freeze out momentum distribution as well as the post freeze out distribution 
are both local thermal equilibrium distributions boosted by the local collective flow velocity on the actual side of the 
freeze out hypersurface, although, the post freeze out distribution need not be a thermal distribution! 

The case of freeze out across a hypersurface with space-like normal shows this clearly because is time-like and 
da M is space-like, thus p^da^ can be both positive and negative. I.e., p^ may point now both in the post- and pre- 
FO directions. Thus the integrand in the above integral (jl|) may change sign in the integration domain, and this 
indicates that part of the distribution contributes to a current going back, into the front while another part is coming 
out of the front. On the pre freeze out side is unrestricted and p^da^ may really have both signs, because we 
may assume that the freeze out front has a certain thicknessu, and due to internal rescatterings inside this front a 
current is fed back to the pre freeze out side to maintain the thermal equilibrium there. 

On the post freeze out side, however, the distribution, fp 0 (x,p ; da must vanish for momentum four-vectors, p^ , 
which point back into the pre FO direction, i.e. do not satisfy the condition, p^da^ > 0.0’Q Thus, this distribution 
cannot be a Jiittner- or other ideal gas distribution]^] 

Nevertheless, the above conservation laws, have to be satisfied, even if the post freeze out distribution is not a 
local thermal distribution. Since, the kinetic definitions of the energy-momentum tensor and conserved current (s) 
are reliably applicable, the conservation laws across a small element of the freeze out front take the form: 


f*Fo( x iP'i T i n i u ''i d(jl ) P^] dap= [ Ng(x) dap , 
P / Js 


( 4 ) 


* Note, that the contravariant normal when becomes space-like, da 11 , should point into the pre-FO direction to satisfy the 
condition, p^da^ > 0, while the covariant normal, dcr M , always points into the post-FO direction! Thus, the direction of 
the contravariant normal d<r M , in the Cooper-Frye formula goes continuously over from pointing to the pre-FO direction to 
pointing to the post-FO direction while the covariant normal of the FO surface stays directed always in the post-FO direction 
when it goes continuously over from time-like to space-like. 
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( 5 ) 


Is (/ fFo( x ,P'’ T , n ’ ua , da ' y ) P fl P v ^j dv»=j Tg v (x) dcTp . 

Here, the matter is characterized by Tq V and on the pre freeze out side of the front. 

The construction of the post freeze out distribution, fp 0 , is a problem in the case of freeze out fronts with space-like 
normal. _ 

For cut Jiittner distribution the conserved currents were evaluated in refJj. Thus, if we know the 5 parameters 
of the pre FO flow and the local freeze out surface from kinetic considerations, then assuming that the post FO 
distribution, fpoiPiX), is a cut Jiittner distribution, we can completely determine the parameters of the post FO 
matter from the conservation laws (|4[J^) . Although, this way we would formally satisfy the conservation laws and 
we would eliminate the particle current pointing back to the pre FO matter, the strange shape of the cut Jiittner 
distribution makes it difficult to accept it as a physical post FO momentum distribution. 


2. Freeze out distribution from kinetic theory 

Following the ideas introduced in ref.l we can calculate the kinetic freeze out distribution based on four-volume 
emission models. The proposed model, on the other hand, requires extended numerical calculation, so here we intend 
to study some overly simplified models, which might give us some hints about the expected shape of post freeze out 
distributions. 

The freeze-out will turn out to be an exponential process, and after about three mean free pathes the amount of 
interacting matter reduces to 5 per cent. Thus, the sharp freeze out layer turns out to be an over-idealization of 
kinetic freeze out in heavy ion reactions, while it is applicable on more macroscopic scales like in astrophysics.^] 

Let us first demonstrate the kinetic model for a drastically oversimplified situation of a plane FO surface. Let 
us assume an infinitely long tube with its left half (x < 0) filled with nuclear mater and in the right vacuum 
is maintained. We can remove the dividing wall at t = 0, and then the matter will expand into the vacuum. 
By continuously removing particles at the right end of the tube and supplying particles on the left end, we can 
establish a stationary flow in the tube, where the particles will gradually freeze out in an exponential rarefaction 
wave propagating to the left. We can move with this front, so that we describe it from the reference frame of the 
front (RFF). 

In this frame, we have a stationary supply of equilibrated matter from the left, and a stationary rarefaction front 
on the right, x > 0. We can describe the freeze out kinetics on the r.h.s. of the tube assuming that we have two 
components of our momentum distribution, ff ree (x,p) and fi n t{x,p). However, we assume only that at x = 0 ffree 
vanishes exactly and /,; nt is an ideal Jiittner distribution (supplied by the inflow of equilibrated matter), while /,; nt 
gradually disappears and // ree gradually builds up as x tends to infinity. We do not assume a priory that fi n t(x,p) is 
an ideal Jiittner distribution for all x, so we will have different FO results depending on the assumed FO mechanism. 


t On the other hand, if kinetic freeze out coincides with a rapid phase transition, like in the case of rapid deconfinement 
transition of supercooled quark-gluon plasma, the short freeze out hypersurface idealization may still be applicable even for 
heavy ion reactions. It is, however, beyond the scope of this work to study the freeze out dynamics and kinetics in this latter 
case. 
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Fig. 1. The freeze-out distribution, ffree{ x ->V)-> in the rest frame of the freeze out front (RFF) calculated from the model presented 
in sect. The momentum is plotted in units of [T], and T = m is assumed. Contour lines are given at values of 0.9, 0.8, 0.7, ... times 
the maximum of ff ree • Here the center of the Rest Frame of the Gas (RFG) is at rest in RFF, = (1,0, 0,0), however, the Eckart 

and Landau flow velocities of the frozen out matter do not vanish! The figures A, B, C correspond to x = 0.02A, 3A, oo respectively. At 
large distances from the initial point of the freeze out process, x —» oo (C), the distribution becomes a cut Jiittner distribution. The 
earlier stages of the freeze out are, however, characterized by asymmetric distributions, elongated in the freeze out direction, x. This may 
lead to a large-pt enhancement, compared to the usual Jiittner assumption as freeze out distribution used in most previous calculations. 


Let us take first the most simple kinetic model describing the evolution of such a system. Starting from a fully 
equilibrated Jiittner distribution the two components of the momentum distribution develop according to the coupled 
differential equations: 


d x f int (x,p)dx = -0(p M dcrp 


( cos Of 
A 

cos Of, 


fint(x,p)dx, 


d x ff ree (x,p)dx = +0(p M ri<7 M ) p fint(x,p)dx. 


( 6 ) 


Here the interacting component, fi n .u will deviate from the Jiittner shape and the solution will take the form 


fint(x,p) = f Juttner (x = 0,p)exp 


-Q(p^da fJ ,) C0S ^ p x 


( 7 ) 


This solution is depleted in the forward Redirection, particularly along the aeaxis. Inserting it into the second 
differential equation above, leads to the freeze out solution: 


f'free(x,p) — fjuttner(x — 0, Rf) 


1 1 — exp 


-Qfp^da^) 


cos 0^ 
A 



( 8 ) 


At x —> oo this distribution will tend to the cut Juttner distribution introduced in the previous section, (see Figs. 
1, 2, and 3.) The remainder of the original Juttner distribution survives as fi n t , even if x —► oo. In this model the 
particle density does not change with x, barely particles moving faster than the freeze out front (i.e. p^dcr^ > 0) are 
transferred gradually from component /j„ t to component ffree- This is a highly unrealistic model, indicating that 
rescattering and re-thermalization should be taken into account in fi nt . This would allow particle transfer from the 
’’negative momentum part” (i.e. p^da^ < 0) of fi n t to ffree, which is not possible otherwise. 
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Fig. 2. The same as Fig. 1, except here the center of the Rest Frame of the Gas (RFG) is not at rest in RFF, u FFG = ( 7 , 0.5, 0, 0). 

At large distances from the initial point of the freeze out process, x -> 00 (C), the distribution becomes a cut Jiittner distribution, but 

less than half oLthe distribution is cut off! Note that the boosted Jiittner distribution became Lorentz elongated and asymmetric (see 
Fig. 2.10 of refjj.) The earlier stages of the freeze out, here also, are characterized by asymmetric distributions, elongated in the freeze 
out direction, x. 


3. Freeze out distribution with rescattering 


The assumption that the interacting part of the distribution remains the distorted (after some drain) Jiittner dis¬ 
tribution, is of course highly unrealistic. Rescattering within this component will lead to re-thermalization and 
re-equilibration of this component. Thus the re-equilibration and the drain terms are in competition and they 
mutually determine the evolution of the component, fi n t- 

To include the collision terms explicitly into the transport equations, (jo!) leads to a combined set of integro- 
differential equations. We can, however, take advantage of the relaxation time approximation to simplify the de¬ 
scription of the dynamics. 

Then the two components of the momentum distribution develop according to the coupled differential equations: 


f^xfint (•£, p)dx 


cos 0 -* 

= -9(p M do>) ^ p fint (x,p)dx 

+ lfeq(x,p) - fint (ft, p)) ^dx, 


d x ffree{x,P)dx = 


cos 0 -* 

+<d(p fJ, dafj,) — 

^ fint (x,p)dx. 


(9) 


Here, the interacting component of the momentum distribution shows the tendency to approach an equilibrated, 
Jiittner type, distribution with a relaxation length coefficient, X 1 ~ A. Of course due to the energy, momentum 
and conserved particle drain, this distribution, f eq (x,p) is not the same as the initial Jiittner distribution, but its 
parameters, n eq (x), T eq {x) and u^x), change as required by the conservation laws. 
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Fig. 3. The same as Fig. 1, except here the center of the Rest Frame of the Gas (RFG) is not at rest in RFF, = ( 7 , —0.5,0, 0). 

At large distances from the initial point of the freeze out process, x —» 00 (C), the distribution becomes a cut Jiittner distribution, but 
more than half of the distribution is cut off! The earlier stages of the freeze out, here also are characterized by asymmetric distributions, 
but these are not elongated in the freeze out direction, x. 


Conservation Laws. In this case the change of the conserved quantities caused by the particle transfer from 
component int to component free can be obtained in terms of the distribution functions as: 

dN‘ u =-L ( dacos dp/ int {x,p) 

A J po 


and 


dT? v = -L [ —p^p-'Q( p ^da l j,) cos 9pf int (x,p). 

A J po 

If we do not have collision or relaxation terms in our transport equation then the conservation laws are trivially 
satisfied. If, however, collision or relaxation terms are present these contribute, to the change of T^ and N^, and 
this should be considered in the modified distribution function fi n t{x,p). 

b. Immediate re-thermalization limit. As a first approximation to the solution of eqs. (|o[) let us assume that 
X' — > 0, i.e. we have immediate re-thermalization after every step dx. Thus the drain is always happening from a 
component of shape f eq (x,p), with parameters, n(x), T(x) and u FFG (x), and we can assume that fi n t = f e q is of 
spherical Jiittner form at any ^.including both positive and negative momentum parts. Above and henceforth the 
notation is similar to the one inB: h = 8nT 3 e^/ T (2Trh)~ 3 , a = so that ) = ha 2 K 2 (a )/2 is the invariant 

scalar density of the symmetric massless Jiittner gas, b = a/y/l — v 2 , v = dao/d<j x , A = (2 + 2b + b 2 )e~ b , and 

ICn(z, w) = z ~ n t°°d X (x 2 - z 2 ) n ~ 112 e~ x , 

(Zn)i J w 

i.e. KL n (z, z) = K n (z). 

In this case the change of conserved quantities due to particle drain or transfer can be evaluated for an infinitesimal 
dx. We assume that the 3-flow is normal to the freeze out surface, and for simplicity we assume v > 0. In this case 
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the change of the conserved particle currents in the RFF is given by 

dNf = -f 4 ^ [ bK ^ b ) + b ^ 2 - l) 7 2 ( 2 A^(a) - AC,(a, b )) + 

+ jv 2 b 2 (2 Kq(o) - /C 0 (a, b)) + 2v 3 j 3 (b + l)e _b ] , 
dN i = -x W K( 3t;2 “ l)7 3 &(2Ad(«) - ACi(a,6))+ 

+(2 + n 4 7 2 6 2 )(2A' 0 (a) - /C 0 (a, 6))- 
—2A'o(6) + 2u7 2 e _fc [n 2 7 2 (6 + 1) + v 2 b — 1]} , 
and for the change of the energy - momentum tensor in the RFF we obtain that 
dT™ = -f {v 2 1 2 b 2 (3 + v 2 )(2K 2 (a)-lC 2 (a,b))+ 

+(v 2 b 2 — v 2 — l)7&(2Afi(a) — /Ci(a, b))— 

—b 2 (2K 0 (a) - /C 0 (a, b)) + ZK^b) + a 2 K 0 (b)+ 

+v'y 2 e~ b (l+3v 2 )"f 2 A(b) — (2+v 2 b 2 )(b+l) + u 4 (l+ 2 g-) 7 2 6 3 |, 
d T ? x = -f^{^b 2 (2K 2 (a)-)C 2 (a,b))+ 

+vab 2 (2K 1 (a) - ACi(a, 6)) + v 2 ^ 2 b (-a 2 + ±4f^6 2 ) - 6 2 + 

+(f 2 + 3) 7 2 A(&)] e- b }-^dA? , 

d2? x = -XT W 3 + ^ 2 )a 2 (2A^ 2 (a) - AC 2 (a, &))+ 

+u 3 a 3 (2ATi(a) - /Ci(a, 6)) + ^.(3 + i> 2 )6 3 + 

+a 2 (u 4 6 - 1) + (3 c 2 + 1 )A(b)] e~ b } - ^dNf , 

dT™ = -f [-(^ 2 + l)a(2AT 1 (a) - AC, (a, &))- 

—v 2 a 2 (2A' 0 (a)—/C 0 (a, 6)) + ^r(6) - 2u(&+l) e - b ] + ^diVf 

and dT zz = dTf v . Note that in RFF the flow velocity of the re-thermalized component is Uj RFG (x) = 
7<t(*) (1) v { x )i 0, 0)|/{ff, where 7 CT = 1/Vl - v 2 . 

The new parameters of distribution fi nt , after moving to the right by dx can be obtained from dN d and dT^ v . 
The conserved particle density of the re-thermalized spherical Jiittner distribution after a step dx is 

hi(x + dx) = hi(x) + dhi{ x) = yj Ad 1 (a; + dx)Ni F {x + dx) 

where the expressions are invariant scalars. After straightforward calculation the differential equation describing the 
change of the proper particle density is: 

dhi(x) = u? rfg (x) dN itfl (x). (10) 
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Although this covariant equation is valid in any frame, we can calculate it in the RFF, where the values of dN-'s 
were given above. Note again, that the particle drain from fi n t.(x), described by dN^ is constrained to the ” positive 
part” in the momentum space, but after re-thermalization we attribute this to the change in the complete spherical 
Jiittner distribution, fi nt (x + dx). Thus, in order to conserve momentum, we have to obtain a decreased Eckart flow 
velocity after the infinitesimal particle drain. 

For the re-thermalized interacting component Eckart’s flow velocity is the velocity of the RFG, which changes 
with x, so we can actually denote this frame as RFG(a;). For the spherical Jiittner distribution the Landau and 
Eckart flow velocities are the same, E RFG {x) = uf L RFG {x) = u i rfg( x )- Thus we can evaluate the flow velocity 

U i,RFG( X + dx ) : 


U z,RFg( X + dx ) = N i( x + dx )/lJ N i N i^, 
which leads to the following covariant expression 


^ U ^i,E,RFG (*^) 




dN i<v (x) 

hi{x) 


( 11 ) 


where Af v (x) = g^ v — RFG (x) u” RFG (x ), is a projector to the plane orthogonal to uf RFG (x). This equation 
is valid in any reference frame, nevertheless we know the four-vectors on the r.h.s. in the RFF explicitly. Then 
the new flow velocity of the matter evaluated according to Eckart’s definition is u^ E RFG (x + dx) = uf RFG (x) + 

d u i,E,RFG ( X )• 



Fig. 4. The same as Fig. 1, except here the center of the Rest Frame of the Gas (RFG) is not at rest in RFF, = ( 7 , 0.5,0.5, 0). 

At large distances from the initial point of the freeze out process, x —» 00 (C), the distribution becomes a cut Jiittner distribution, which 
is not centralized in p y and less than half of the distribution is cut off! The earlier stages of the freeze out, here also are characterized by 
distributions asymmetric both in the directions p x and p y , and these are also elongated in the direction of the freeze out flow velocity, 

,, v 

U RFG' 

To get the temperature and the change of Landau’s flow velocity, we have to analyze the change of the energy 
momentum tensor. Before the particle drain the energy - momentum tensor at x in the RFG is diagonal, T/^(x) = 
diag(e*, P i: P z , Pi)\ RFG {x) while in the RFF T? v {x) = [(e* + P t ) u^ RFG u^ RFG {x) -P^g\ RFF ( X ). Adding the drain 
terms, dT^ v (x), to this arising from the freeze out while we move to the right by dx, yields (x + dx) which will 
not be diagonal in the RFG(a:) and the pressure part will not be isotropic. We can Lorentz transform this to 




another frame which diagonalizes Tj 1 (x + dx). This means to find the Landau flow velocity of the new system, 
u iL rfg( x + dx ) original RFG(a;). After a straightforward diagonalization, a somewhat tricky algebra and 

neglecting second and higher order terms we arrive at the covariant expression^ 


du i,L,RFG ( X ) 


AHx) dT iF<? <,rfg ( x ) 

&i + Pi 


( 12 ) 


Although, for the spherical Jiittner distribution the Landau and Eckart flow velocities are the same, the change of 
this flow velocity when calculated from the baryon current and from the energy current are different 


du i,E ,rfg ( x ) ^ du i,L ,rfg ( x ) ■ 

This is a clear consequence of the asymmetry caused by the freeze out process as we pointed out already at the 
discussion of the properties of the cut Jiittner distribution. Unfortunately, this also illustrates the weakness of our 
assumption on the complete re-thermalization to a spherical Jiittner distribution, because we cannot choose the 
correct velocity change: If we choose du ^ as the new velocity of the (spherical Jiittner distribution) fi n t{x + dx), 
then we violate the momentum conservation in our model, on the other hand if we choose du^, then we violate the 
baryon current conservation! Thus a spherical (or even elliptic) distribution cannot be fitted to the freeze out drain, 
and we would have to use an ansatz, which has (in addition) an asymmetry in the x direction (i.e., an egg shape), 
for the distribution fi nt . 

Being aware of this weakness of the model, we nevertheless, maintain the assumption of spherical Jiittner shape 
for fi nt for the sake of simplicity. We can choose the flow velocity change then according to the physical problem. 
For example for the freeze out of baryon free plasma this problem does not occur, and we have to choose du 

The last item is to determine the change of the temperature parameter of /j„ t . From the relation e = u IJ T tlv u l , 
we readily obtain the expression for the change of energy density 


dei(x) = u^i,RFG{x) d.T? v (i r) u Vti ^ RFG {x ), (13) 

and from the relation between the energy density and the temperature (see Chapter 3 in refB), we can obtain the 
new temperature at x + dx. Fixing these parameters we fully determined the spherical Jiittner approximation for 
fint ■ With this ansatz the pressure asymmetry and pressure balance cannot be realized, thus our model will be only 
a rather approximate description of the freeze out process. 

Nevertheless, we can draw some preliminary conclusions about the development of the kinetic distribution during 
freeze-out. 


4. Conclusions 

We turned to the problem of estimating the freeze out distribution. Obviously the real freeze out distribution depends 
strongly on the details of the freeze out (and hadronization) dynamics. In heavy ion reactions, the curvature of the 
freeze out surface and the conditions varying in time do affect the freeze out distribution, nevertheless, as a first 
step, we assumed that the process is stationary and the curvature of the front is negligible. These approximations 
are extreme, but still enable us to draw some preliminary conclusions. 


* Let the energy-momentum tensor of a system be T pv . The energy and momentum flow is characterized by the Landau flow 
velocity, a unit four vector, u p . We are looking for a relationship between the infinitesimal change of the flow velocity du p 
and the cqsresponding shift in the energy-momentum tensor We introduce the projector, A MI/ = g — ii M iA, with the 

properties Q A^Uu = 0 and du p = A p du„ since Updu p = 0 . The Landau flow velocity is parallel to the flow of the momentum. 
Thus Up = Const, x T p u v , therefore AppT^Uv = 0. We differentiate the above equation and take into consideration the 
identities e = and A p T ptT A" = —PA 111 ', where e and P are the energy density and pressure of the dissipationless, 

fully equilibrated fluid. Then using the properties of A pv we get du p {e + P) + u p du ll T ,J '''u l , = A P pdT pv u v . Since the flow 
velocity and the momentum flow are parallel the second term on the l.h.s vanishes. Thus the equation describing the change 
of Landau’s flow velocity becomes du p = A PPb dT pv u v /[e + P). 
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Following the lines and ideas presented in rebel, the first simple kinetic freeze out model reproduces the cut Jiittner 
distribution as the limiting distribution, ff ree , after complete freeze out at large distances. However, the model 
at the same time leads to unrealistic consequences, namely that the interacting part of the distribution, /, nt also 
survives fully, as the other part of the Jiittner distribution. Thus having both components at the end in this model, 
the physical freeze out is actually not realized. This turns out to be a consequence of the fact that the effect of 
rescattering and thermalization in the interacting part of the distribution was ignored. 

In an improved but still rather approximate kinetic freeze out model which takes rescatterings into account, the 
interacting component is assumed to be instantly re-thermalized taking a spherical Jiittner shape at each time step 
with changing parameters. The model leads to a set of coupled differential equations (|lO|Jl]]Jl^ , |l3| ). Equations ([□]) 
and © can be used in some combined form, or one of them can be selected which fits the physical situation the 
best. Then the three parameters of the interacting component, fi n ti can be obtained in each time step analytically 
(considering IC n (x,y) an analytic function). 

Now the density of the interacting component will gradually decrease and disappear according to eq. (10), the flow 
velocity will also decrease in both cases, a or because only forward going particles freeze out, and the energy 
density will decrease also according to eq. ©■ Thus, the initial contribution to ffree at small x will resemble the 
distribution shown in Fig. 2A, then as x increases and the velocity decreases it will become to similar to Fig. IB, 
while at the final stages it will approach Fig. 3C. As a consequence the integrated distribution will not resemble a 
cut Jiittner distribution. 

Thus the arising post freeze out distribution, ff ree will be a superposition of cut Jiittner type of components, 
from a series of gradually slowing down Jiittner distributions. This will lead to a comet shaped final momentum 
distribution, with a more dominant leading head and a tail. In these rough models a large fraction (~ 95%) of the 
matter is frozen out by x = 3A, thus the distribution ffree at this distance can be considered as a first estimate 
of the post freeze out distribution. One should also keep in mind that the models presented here do not have 
realistic behavior in the limit x —> oo, due to their one dimensional character. Nevertheless, this improved model 
with rescattering enables complete freeze out (unlike the simpler model in sect. where only the originally forward 
moving particles freeze out even at large distances). 

In case of rapid lradronization of QGP and simultaneous freeze out, the idealization of a freeze out hypersurface 
may be justified, however, an accurate determination of the post freeze out hadron momentum distribution would 
require a nontrivial dynamical calculation. 
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